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ABSTRACT 

We examine the properties and evolution of a simulated polar disc galaxy. This galaxy 
is comprised of two orthogonal discs, one of which contains old stars (old stellar disc), 
and the other, containing both younger stars and the cold gas (polar disc) of the 
galaxy. By exploring the shape of the inner region of the dark matter halo, we are able 
to confirm that the halo shape is a oblate ellipsoid flattened in the direction of the 
polar disc. We also note that there is a twist in the shape profile, where the innermost 
3 kpc of the halo flattens in the direction perpendicular to the old disc, and then 
aligns with the polar disc out until the virial radius. This result is then compared to 
the halo shape inferred from the circular velocities of the two discs. We also use the 
temporal information of the simulation to track the system's evolution, and identify 
the processes which give rise to this unusual galaxy type. We confirm the proposal 
that the polar disc galaxy is the result of the last major merger, where the angular 
moment of the interaction is orthogonal to the angle of the infalling gas. This merger is 
followed by the resumption of coherent gas infall. We emphasis that the disc is rapidly 
restored after the major merger and that after this event the galaxy begins to tilt. 
A significant proportion of the infalling gas comes from filaments. This infalling gas 
from the filament gives the gas its angular momentum, and, in the case of the polar 
disc galaxy, the direction of the gas filament does not change before or after the last 
major merger. 
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1 INTRODUCTION 

Classical disc galaxies consist of flattened distributions of 
stars, gas, and dust, co-rotating with the angular momen- 
tum vectors more or less alignedQ polar disc galaxiefl are 
significantly more rare, possessing two orthogonal discs. An 
original survey of galaxies con sisting of orthogonal s truc- 
tures of stars was carried out bv lWhitmore et all (| 1 QQdl ') who 
identified a sample of polar ring galaxies. The 'primary' discs 
in such systems are typically classified as lenticulars (SO), or 
Ellipticals, with a second, younger, structure aligned with 
angular momentum axis of the primary disc. Fewer than 1% 
of SO galaxies contain a polar structure. 



1 In the presence of warps thoug h, there may be so me misalign- 
ment driven by external torques (Roskar ct al. 2010). 
Originally called 'multi-spin galaxies' bv lRubinl lll994h . 



Subsequent studies of NGC4650A revealed a subclass of 
structures where the orthogonal distribution of stars was not 
a ring but a second fully-formed continuous disc structure, 
Jlodice et al.ll2002l: [Gallagher et al1l2002l; ISwaters fc Rubinl 
120031 ; iMaccio et al.1 l200fj ; iBrook et all boOgP . Systematic 
studies suggest that their stellar populations, light distri- 
bution, and gas-phase characteristics, are not dissimilar 
from those of cl assical discs - e.g., exponential light profiles 
jSchweizer et al.1 1 19831 ) and HI mass-to-optical luminosities 
typi cal of mass to B-band luminosity typical of late- type spi- 
rals (lHuchtmeier||l997l : lArnaboldi et al.lll997l ; ISparke fc Coxl 
2000 : ICox et alj|2006l ) . 



Thus there exist polar ring galaxies as well as polar disc 
galaxies. We will discuss the latter type in the following pa- 
per. On of the important uses of polar disc galaxies is a 
direct result of having orthogonal discs at the bottom of 
a halo potential well. One stellar disc historically allowed 
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astrono mers to calculate the e xistence of dark matter ha- 
los, (see Begeman et al] |l99ll ) for a detailed discussion). 
Ilodice et alT l|2003) note that the existence of two orthog- 
onal discs means that the shape of the inner region of the 
halo can be recovered from observations. 

There are several theories proposed to explain the origin 
of polar ring and polar disc galaxies: 



(i) Mergers : As e nvisio ned by iBekkil l| 19971 . ll998T l and 
IBournaud fc Combes! l|2003l ). polar rings can eventuate for 
very specific collision configurations where the collision is 
'head-on' and the initial angular momentum low. 

(ii) Tidal Accretion : Proposed by ISchweizer et al.l 

and simulated by iReshetnikov fc Sotnikoval l|l997h ~ 
this scenario involves the capture of a gas-rich secondary by 
the (future) polar disc host halo. The interaction-induced 
polar ring is formed more rea dily than in the aforemen- 
tioned merger-induc ed scenario (IBournaud fc Combe s 2003; 
ICombes et alll2006h . 

(iii) Infall from Filaments : If gas falls into a galaxy 
along cosmic filaments that are inclined to the stellar disc, a 
polar d isc may form jMaccio et al.l 120061') fc llBrook et al. 
I200SI) . ICombes et ail (|2006t ) and IBournaud fc Combes 
(2003) suggest such cosmological infall is a viable alter- 
native to tidal ac c retion for the NGC 4650 A system e.g. 
l|lodice et al.ll2006l ; ISpavone et aT]|2010l . l201ll ). 

(iv) Resonance : This approach assumes that the polar 
disc is formed via resonant coupling to a triaxial halo poten- 
tial |Tremaine fc YulbOOOT ). In this picture, a disc is taken to 
lie in the plane of symmetry of the triaxial halo and the halo 
tumbles with respect to the disc. As the tumbling slows, the 
stars in the disc can get trapped in resonance with the halo. 
A stellar orbit inclined to the disc precesses at a slow ret- 
rograde rate and the star can be pushed into a polar orbit. 
This method has been invoked to best explain equal mass 
stellar disc systems. 



A seminal rev iew of the field is provided by 
ICombes et alj (|2006t ). to which all interested readers should 
refer for a rich surve y of the field. The fil amentary infall sce- 
nario (iii, above) of iMaccio et al.l (|2006h has gathered mo- 
mentum over the past few years as being the favoured origin 
of polar disc systems. This can be traced, in part, to the fact 
that the alternatives do not account fully for the range of 
observations - in particular, those which show that the polar 
disc is often of high mass, possesses extended and continu- 
ous structure, and lacks obvious evidence of starbursts that 
might be associated with polar structure formation. Scenario 
(i) can form a polar disc/ring gala xies when the relative i ni- 
tial velocities are low (e.g. IBournaud fc Combes] (|2003h ) 
and by (ii) in loose groups of galaxies such as UGC9796 
(jSpavone et alJ[201lT ). Polar disc/ring galaxies are found in 
low density e nvironments as the polar structure is destroyed 
by mergers (|Macci6 et al] 120061 ) . This suggests that each 
type of formation may exist in the universe. 

Cl osely related to | Macci6 et al] ((2006) scenario (iii, 
above). [Brook et all (|200sh suggest that the polar disc forms 
as a result of the last major merger changing the angular 
momentum of the stars and gas of the stellar disc, and then 
subsequent gas continues to fall in along the old trajectory 



and builds up the polar disc[f] The basic observable charac- 
teristics of their simulated polar disc system were presented, 
including star formation rates, circular velocity profiles, and 
structural parameters of the discs, along with the aforemen- 
tioned putative origin scenario. The exact mechanism by 
which the dark halo aligns to the orthogonal discs and how 
the new infalling gas be comes misaligned to the orthogonal 
discs was not studied bv lBrook et al] (|2008T ). 

In what follows, we rectify this by examining the un- 
der lying physics gover ning the misalignment of the discs in 
the lBrook et al] (|2008T ) simulation, tracing its temporal evo- 
lution and association with large-scale structure, in order to 
put their origin scenario to the test, quantitatively. The sim- 
ulation and its basic properties are reviewed in § 2 and § 3, 
while the time evolution of the host halo's shape is detailed 
in § 4. The metallicity and age gradients of the discs are 
then confronted with recent observational work (§ 5), while 
the time evolution of the dark halo, gas, and stellar angular 
momenta is derived in § 6. Our conclusions and suggestions 
for future directions are made in 8 7. 



2 CODE 

The polar disc system under examination here was realised 
with the gravitational N-body + smoothed particle h ydrody- 
namics (SPH) code Gasoline (IWadslev et alj|2004l) . Grav- 
itational forces are computed using a iBarnes fc Hut] (|l986l ) 
tree-based algorithm; the code is spatially and temporally 
adaptive and employs periodic boundary conditions. SPH 
is used to discretise the astrophysical gas and accounts 
for radiative cooling, ionisation, and excitation of the gas. 
Th e ultraviolet background radiation field is included, af- 
ter |Haardt3kJMadau| (Il996l ). The internal energy of the gas 
is calculated to conserve energy better than entropy and 
shocks are accounted for via artificial viscosity. 

The st ar formation re c ipe em ployed within Gas oline is 
detail ed bv lStinson et al] l|2006h and includes (i) a iKroupal 
l|200ll ) initial mass function (IMF), (ii) the effects of both 
Type la and Type II supernovae (neglecting the instanta- 
neous recycling approximation), and (iii) a blast- wave feed- 
back scheme with which to couple supernova energy to the 
surrounding interstellar medium. 

F or this simulation, first introduced by iBrook et al] 
(2008), volume re-normalisation was employed to ensure 
high spatial and force resolution in the area of interest, 
while self-consistently tracking the gravitational torques im- 
parted by large-scale cosmological structure. The cosmo- 
logical framework was a concordant ACDM, with fio=0.3, 
f2A=0.7, o"8=0.81, h=0.7, where Qq is the matter density, 
Qa is the fraction of the energy density of the Universe due 
to the cosmological constant, ag is the rms density fluctu- 
ation in 8 Mpc sub-samples, and h the Hubble Constant. 
The initial power spectrum was calculated using C mbfast 
|Zaldarriaga et alj[l99sl ; ISeliak fc Zaldarriagalll996h . 

An initial, coarse, dark matter only simulation was run 
to redshift 2=0 in a 28.5/i _1 Mpc box; this is large enough 
to provide realistic tidal torques with sufficient resolution 



3 Alternatively, the gas falling along filaments could change di- 
rection. 
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Figure 1. Schematic density distributions of the intermediate- 
age stars associated with the primary stellar disk (colours; vertical 
in this plane) and cold gas associated with the secondary polar 
disk (contours; horizontal in this plane). 



for halo identification. A host halo was selec ted using AHF 
l|Gill et al.ll2004lKnollmann fc Knebe!l2009l ), and the parti- 
cles therein traced back to their respective locations in the 
initial conditions. The volume is then re-centred on this re- 
gion and said region 're-populated' with lower mass dark 
matter particles in concentric shells of increasing resolution. 
Each 'generation' of dark matter particles is 1/8 the mass of 
the previous generation. The highest resolution dark matter 
particles have a mass 9.5 x 10 4 Mq. Only within the high 
resolution region are gas particles now placed, and tracked 
forward in time in a re-simulation of the volume, now in- 
cluding the effects of star formation and feedback. These gas 
particles each have a mass 1.6 x 10 4 Mq; star particles which 
form from these gas particles have masses of 3.3 x 10 3 M . 
The force resolution within the inner, high-resolution, re- 
gion is 150 pc. At redshift z=0, there are 3.5 x 10 6 baryonic 
particles within the inner 20 kpc. 



3 BASIC CHARACTERISTICS OF THE 
SIMULATED SYSTEM 

As noted earlier, the simulation discussed here was first in- 
troduced by Brook et al. (2008), but a detailed examina- 
tion of its dark matter halo characteristics, relationship to 
large scale structure, stellar and gas-phase metallicity dis- 
tributions, and their relationship to the system's proposed 
formation, was reserved for future work. We address these 
points in what follows. A schematic of the system is shown 
in Fig. [T]where the older, 'primary', stellar disc is shown in 
colour (vertical in this figure) and the cold gas associated 
with the younger, 'secondary', polar disc is represented by 
the red contours (horizontal in this figure). 

We can identify the polar and stellar discs cleanly via 
the ages of their associated star particles : recently-formed 
stars (t >9 Gyrs) define the 'polar disc', while intermediate- 
age stars (4< t <6 Gyrs) define the 'stellar disc'. This simple 
delineation is used to define the disc orientations for our 
subsequent analyses. Stars formed early (t <4 Gyrs) are 
associated with the halo of the galaxy and can be traced 



to the stars which formed before or during the last major 
merger. 

The system as a whole has a dynamical (stellar) mass 
of 1.6 x 10 11 M (1.0 x 10 10 M Q ) within its virial radius 
(~130 kpc) at z=0.17. The stellar mass of the old stellar 
disc is 1.9 x 10 9 Mq while the polar disc has a mass of 9.8 x 
10 8 Mq - i.e., the polar disc has half the stellar mass of the 
stellar disc. The gas mass in the polar disc is approximately 
the same as the stellar mass at z=0.17, so the baryonic mass 
in the two discs is approximately the same. Both discs are 
less massive than the very old spheroid, (with stars formed 
at less than 4 Gyr which has a mass of 6.2 x 10 9 Mq). 
Of the three polar ring galaxies studied in Spavone et al. 
(2010, 2011), two show a polar disc more massive than the 
host galaxy, while UGC 7576 has a polar ring 2.7 times less 
massive than the host galaxy. Our polar disc is 3.6 times less 
massive than the host, a valu e similar to that found for UGC 
9796 bv lSpavone et all (|201ll ) The simil arity in the masses o f 
the two discs differs from NGC 4650A (|Spavone et al.ll201lD 
These mass values are highly dependant on the age cuts 
used t o define the disc. The original values from lBrook et all 
(2008) give significantly different mass values for the two 
discs and so the relative masses of the two components is 
mostly dependant on how one chooses to define it, rather 
than an intrinsic pro perty. For example , if we decompose 
using the method of lAbadi et al.l (|2003l n then the stellar 
masses of each disc are the same, so that the total baryonic 
mass of the polar disc is twice that of the stellar disc. 

The stellar disc is essentially gas-free by redshift z~0.2, 
but prior to this, there is an extended period of time dur- 
ing which there is gas in both discs. The scale-lengths of 
the stellar and polar discs are 1.0 and 3.4 kpc, respec- 
tively, with both showing disc truncations near ~3 disc scale- 
lengths. The two discs are nearly orthogonal, with an incli- 
nation of ~84° between the two. The circular velocity at 
2.2 scale-lengths fo r both discs is ~110 km/s, based on the 
iBrook et al. (12008 1 ) ana lysis. The polar disc system reported 
bv lMaccio et al.l (|2006l ) possesses similar masses and radial 
extents for its two discs. We show the mass accretion his- 
tory of the galaxy in Fig. [2] We can, in some sense, think 
of this history in terms of three epochs of 'star formation' 
: a rapid initial phase (t<4 Gyrs), an intermediate phase 
(4<t<8 Gyrs), and a more gentle, later, phase (t>8 Gyrs). 



3.1 History of Formation 

The primary disc first forms at redshift z=3.3 (age=2 Gyr) 
and the last major merger (with a mass ratio near 1:1) 
encountered at z=2.2 (age=3 Gyr). The interaction be- 
tween the merging galaxies lasts for ~1 Gyr and signif- 
icantly disrupts the primary disc, although the disc re- 
establishes itself rapidly due to its gas-rich nature. There 
is a delay between the last major merger and the time 
at which the polar disc becomes noticeable (at redshift 



4 This galaxy is decomposed according to the ratio of e = 
Jz(E)/ J c i rc (E), where Jz(E) is the z component of the angular 
momentum of the galaxy and J c irc(E) is the angular momentum 
of the circular orbit with the same energy E. Disc stars are chosen 
via a cut in e close to 1. 
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Figure 2. The upper panel shows the mass accretion of the bary- 
onic components, with galaxy mass against the age of the uni- 
verse. The bottom panel shows the d ark matter accretion h istory, 
where the dotted line is the best fit (Wcchslcr ct al. 2002) func- 
tion. 
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Figure 3. Time evolution of the star formation rate of the sim- 
ulated system. The solid black line shows the star formation rate 
inferred from stars presently situated within a slit aligned with 
the old stellar disc and excluding the first kpc; the solid red line 
corresponds to the star formation rate inferred from a slit aligned 
with the polar disc; the dashed line corresponds to the global star 
formation rate of the syst em scaled down by a factor of ten. This 
is an expanded version of lBrook et al. figure 2. 



0=0.8, age=6.8 Gyr). The polar disc then establishes it- 
self, evolving passively and stably f o r the last ~4 Gyrs of 
the simulation. In the I Brook et al.l |200ct ) galaxy studied 
here, the polar disc component of the galaxy initially forms 
as a ring but over the following Gyrs the centre is filled 
by cold gas to resemble a disc. For a video o f the polar 
disc formation, first presented in iBrook et al l (|2008l ). see 
http:/ /www.youtube.com/watch?v=c-H3WzaewdY Only 
in the final ~0.2 Gyr does the polar disc become unsta- 
ble due to coupling between the two discs, and the polar 
disc structure vanishes when the two align. There is some 
evidence of a bar growing in the polar disc right up until 
the disc is disrupted, and a bar is evident in the old stellar 
disc as well. The surprising longevi ty of simulated polar disc 
structures w a s high lighted by both lMaccio et all (|2006l ) and 
IBrook et al.l (|2008l ). The former note that only gas which 
falls into the polar disc at ~90° to the massive stellar disc 
can form a stable secondary disc; at smaller angles of in- 
cidence, the infalling gas was found to fall rapidly to the 
already existing disc. 

3.2 Star Formation History 

The polar disc system experiences a fairly constant, intense, 
phase of star formation between 1—4 Gyrs after the sim- 
ulation commenced (Fig [3}. Subsequent to this, the global 
star formation decreases rapidly for the final ~8 Gyrs of the 
simulation (dashed black curve). We can infer the star for- 
mation histories along each disk by placing 'slits' along each 
(with the knowledge that by doing so means some 'dupli- 
cation' of star formation histories in the region of overlap 
between the slits). Although we use the two discs to define 



where we place the 'slits' we include all stars along the 'slit' 
in our determination of the SFR. The result of this suggests 
that the primary stellar disc forms its stars preferentially at 
early times, prior to the last major merger, although resid- 
ual gas does remain to form stars at later times (solid black 
curve) . Over the final ~5 Gyrs of the simulation, stars start 
to form in the polar disc (solid red curve - see also Fig 2 
of Brook et all (|2008l )1. increasing steadily, while the star 
formation in the primary stellar disk has essentially halted 
over the final ~4 Gyrs. During this era, the star formation 
in the polar disc is very smooth and continuous, with no 
starburst features indicative of merging events. Star forma- 
tion is ongoing and increasing, up to z=0.17 but there is no 
sign of bursting, with a considerable fraction of stars being 
produced recently, so we expect young stars in the polar 
disc, as in NCG 4650A ^Gallagher et al.ll2002r i. although no 
spiral can be seen such as in NGC 4650A. Star formation 
is uniform throughout the disc, so there are no young star 
cluster regions, merely global star formation without clump- 
ing. This does not, however, m ean that the SFR is s imilar to 
UGC 7576 or UGC9796 which lSpavone et al l (|201lD fit with 
an exponentially decaying SFR. Indeed, the SFR in the po- 
lar disc is increasing with time. The residual star formation 
not accounted for between the dashed curve and the sum 
of the two solid curves is that associated with the spheroid 
components of the system. In Fig the dotted line (rep- 
resenting 10% of the global star formation rate) shows that 
the star formation within the discs (the 'in situ' component) 
is a small fraction of the global star formation at early times. 

As noted in § 3, the system was decomposed using a 
stellar age cut : stars formed after t>9 Gyrs are considered 
as part of the polar disc, and those formed before that are 
assigned to the stellar disc. 
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4 THE HALO 

One of the original reasons for postulating the existence 
of dark matter was the circular velocity profiles of galaxy 
discs. With polar discs, we have circular velocities of two 
orthogonal components of the dark matter density profile 
and so these objects provide an ideal method for study- 
ing the shape of the d a rk m a tter halo, as highlighted 

I) ICasertano fc van Gorkonj 
Sackett et al.l [| 19941 ) and 



by Combes fc Arnaboldil 
(1991), ISchweizer et al, 

llodice et alj |2006l ); indeed, making use of the Tully-F isher 
relation, as applied to polar discs, (jlodice et al.ll2003h sug- 
gested that dark matter haloes are oblate and flattened 
along the polar disc because the HI gas polar disc was found 
to have a faster circular velocity than the stellar disc. 

In calculating th e shape of the dark m atter halo directly 
we use the method of iKnebe etafl (|201Ch . The shape of the 
halo is calculated using spherical shells; at each radius, Ri 
we calculate the inertial tensor using all particles within that 
radius. We use the standard inertial tensor, as opposed to 
the reduced tensor, because this technique strongly biases 
the profile in favour of the central region. 

The halo shape is calculated, according to the inertial 
tensor, as 



(i) 



where Iij is the inertial tensor, n is the particle id, i and j 
are the x,y,z components of the particle position vector. The 
axial ratios are (a, b, c) = (\/Ar, <*/\b, y/X^) where a > b > c 
and \ a ,b,c are the eigenvalues of the tensor. The directions 
of the axes are the correspon ding eigenvectors, as outlined 
bv lBailin fc Steinmetzl (120051 ). 

This method easily identifies the direction of the prin- 
ciple axes of the best fit ellipse, but introduces a bias in the 
halo shape. To correct this to match halo sphericity based 
on potential energy contours we must adjust the sphericity 
parameter, S = c/a by, 



>-3tr tie — 'J. 



VS 

sphere ' 



(2) 



where the sphericity is the ratio between the major and mi- 
no r axes of the halo. This re lation is based upon the study 
bv lBailin fc Steinmetzl (|2005l ) 

We will c ompa re our inferred halo shape with that of 
iMaccio et al.l ([2006), in addition to the work of lBailin et al.l 
j|2005l) . who studied the alignment of traditional galaxy discs 
with their dark matter haloes and found that in the inner 
region there was a tight alignment which becomes increas- 
ingly divergent with distance. Thus, we have the amount of 
mass within a radius R along two orthogonal directions, Mi 
and M 2 . 

In order to calculate the axial ratios we imagine two 
spheres of radii 'a' and 'b', containing the same mass. The 
ratio of the volumes of these two spheres is the same as 
the ratio M2/M1. Therefore the axial ratios of the halo are 
(M 2 /Mi) 1/3 . 



4.1 Radial Profile 

Using the inertial tensor approach outlined above, in Fig 4 
we show the sphericity profile of the system's dark halo (up- 
per panel) and the direction of the minor axis of the in- 
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Figure 4. Upper panel: The sphericity profile of the dark mat- 
ter halo as a function of galactocentric distance. Bottom panel: 
The direction of the minor axis of the inertial tensor (black: z- 
component; red: x-component; cyan: y-component). 



ertial tensor (bottom panel) as a function of galactocen- 
tric distance. We find the system is embedded within a 
nearly spherical (S=0.9) halo near the edge of the stellar 
disk (~5 kpc), but which declines in sphericity both within 
the inner kpc and beyond ~10 kpc. This reaches a minimum 
sphericity at 54 kpc or approximately half the virial radius 
and increases again to a final value of 0.67 at R„i r . The dark 
halo is prolate at the edge of the stellar disc, distinctly oblate 
at about 15 kpc and then becomes increasingly triaxial. No 
obvious 'imprint' of the polar disc is seen in the sphericity 
profile. 

This sp h ericity trend is similar to that of the 
M accio et all (|2006l ) galaxy, which has an inner region of 
greater sphericity than the outer region. This will be dis- 
cussed further, later in this section. 

In the bottom panel, we can see that the direction of 
the halo's minor axis in the innermost region is orientated 
90° to the minor axis in the outer region. The inner region is 
flattened against the direction of rotation of the stellar disc 
and the polar disc is aligned to the axis of rotation of the 
stellar disc. The polar disc has been aligned to the z-axis of 
the volume and the stellar disc is approximately aligned to 
the y-axis. Thus, it is clear that the minor axis of the halo 
in the innermost region is aligned to the stellar disc and in 
the outer region it is aligned to the polar disc. 

For non-polar disc galaxies the minor axis of the dark 
matter halo is found to a lign closely with th e angular mo- 
mentum of the disc (e.g., iBailin et al.1 (|2005l ). Fig 5 shows 
an example o f such a non-polar disc galaxy (gl5784 from 
IStinson et ail l|2010l) ). This galaxy shares some similarities 
with our polar disc galaxy in that it has not experienced a 
merger or starburst for the past ~4 Gyrs. From the bottom 
panel, one can see that there is no sign of the obvious twist 
in the minor axis profile (cf. Fig 4, for the polar disc system 
here), consistent with the conclusions of Bailin et al, thus 
demonstrating a significant difference in the properties of 

the two dark haloes. 

Another feature of the dark haloes of the IBailin et al.l 

(2005) sample of eight galaxies is that although the disc 
angular momentum and dark halo minor axis are tightly 
aligned in the centre, they diverge with increasing radius. 
The dark halo of the polar disc galaxy seems to remain 
closely aligned at large radii outside the twisted region. Such 
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Figure 5. Top panel: The sphericity profile of the dark matter 
halo as a function of galactocentric distance, for a simulated 'clas- 
sical' disc galaxy. Bottom panel: The direction of the minor axis of 
the inertial tensor (black: ^-component; red: ^-component; cyan: 
y-component). 



behaviour is seen in only one galaxy of the eight in the Bailin 
et al. sample. This may be a coincidence, but this strong be- 
haviour in an object that is itself uncommon may hold addi- 
tional clues as to the physics by which the polar disc formed. 
A larger sample size is required to exclude the possibility of 

cha nce. 

llodice et al.1 (|2003l ) note that polar disc galaxies have 
an offset from the Tully-Fisher relation relative to other spi- 
ral galaxies, showing faster rotation for a given magnitude. 
This is indicative of flattening of the potential towards the 
plane of the polar disc, as we have demonstrate d above. Our 
polar disc galaxy (Mi=— 20.3; B— 1=1.7, from 
20081 ) Table 1. This is for the whole galaxy as in 



Brook et al 



Iodice et al 



2003)), w hen compared with t he Tully-Fisher relation of SO 
galaxies (|Williams et al.l 120101 ) lies above t he best-fit line, 
but n ot significantly so. Considering that Wi lliams et al.l 
(2010) note that the Tully-Fisher relation for SOs is lower 
than that of spirals, the polar disc galaxy should actually 
be indistinguishable from these standard populations, as we 
find when we co mpare our stel l ar ma ss an d circular velocit y 
to the results o f iPizagno et ail |2005l ) and lMcGaughl (|2005h . 
among others. Thi s is not unexpecte d as only 7 of the 16 
galaxy sampled bv llodice et al.l l|200St ) are exceptionally dis- 
placed from the general Tully-Fisher profile in their paper. A 
larger sample of simulated polar discs would be required to 
develop any further conclusions in this area, together with 
a lar ger sample of K-ban d observations of PRCs. 

iMaccio et al.l ((2006) studied the shape of a simulated 
polar disc galaxy and calculated the sphericity of the galaxy 
halo as we have done here. There are some differences. Their 
shape profile peaks at a sphericity of approximately 0.8, at 
a radius of 9.8 kpc. Our polar disc peaks at a sphericity of 
0.93 at only 3 kpc, corresponding to the edge of the stellar 
disc. At half the vi rial radius our halo i s still more spherical 
than was found bv lMaccio et al.l (|2006t ) with a sphericity of 
0.65 compared to their 0.5. Our values are fairly stable with 
redshift during the lifetime of the polar disc. This galaxy 
has a sphericity peak at approximately twice the inner disc, 
unlike in our galaxy where the peak an d the edge of the stel - 
lar disc coincide. The sphericity peak in lMaccid et al.l (|2006l ) 
is 0.8 compared to our 0.9 and at half the virial radius the 
sphericity is 0.5, again, less spherical than for our polar disc. 



Without further comparison in the shape finding methods, 
and a larger statistical sample further conclusions cannot be 
drawn. 



4.2 Evolution 

We analysed both the classical disc galaxy from Fig. 14.11 
and the polar disc galaxy towards the ends of their re- 
spective quiescent phases (~4 Gyrs after the last significant 
merger/burst of star formation, or post-establishment of the 
polar disc, respectively). The direction of the minor axis of 
the polar disc's dark halo is strongly aligned with the same 
plane (the z-axis of the polar disc angular momentum) for 
much of the simulation (Fig. 14. 2j) . This is not the case for 
the classical disc, where the minor axis changes significantly 
from redshift ,z=0.3 (age=10.3 Gyr) to 2=0. Thus, there ap- 
pears to be a strong, long term alignment in the shape and 
direction of the polar disc halo. This strong alignment may 
be important in giving rise to the polar disc, but we are 
hindered by our sample size of one in coming to a strong 
conclusion on this point. 

The inner region of the polar system shows a long-term 
'twist' that evolves towards orthogonality. The classical disc, 
after the last starburst, sees an evolution towards alignment 
throughout the halo, such that no twist is evident. We sug- 
gest that the initial misalignment is due to the merger of in- 
falling sub-haloes temporarily re-orienting the disc and the 
alignment thereafter is traced to the later quiescent phase. 
This raises the interesting possibility that alignment of the 
galaxy and da rk halo is due to th e torques on the galaxy by 
the dark halo. iBailin et alj (|2005t ) discounted this, although 
they did find that the disc is always aligned with the inner 
halo, which we find to not be the case, at least for part of the 
lives of both discs. For our polar disc system, the orientation 
of the dark matter halo in the outer region is set early and 
does not significantly change over the course of the simula- 
tion. The classical disc shows considerable variation in the 
direction of the halo minor axis from 0° to 90° over time. 



4.3 Shape from observations 

In order to derive the shape of the polar disc halo from obser- 
vational data alone we can take measurements of the circular 
velocity of each disc. This is best done from the edge-on view 
because it becomes difficult to i dentify these galaxies sim- 
ply because of orientation effects |Whitmore et al.1ll990h . By 
measuring the circular velocity of the two discs edge on we 
can gain insight into the shape of the dark halo. 

First, we derive the line-of-sight velocity profiles of the 
two discs by, again, lying down 'slits' along both. These pro- 
files can then be used to infer the enclosed mass as a func- 
tion of galactocentric radius. If the inferred mass along one 
disc is greater than that inferred along the orthogonal slit, 
the former is likely the direction of the major axis. We use 
the aforementioned redshift 2=0.17 snapshot (the point at 
which the two discs are near the end of their ~4 Gyr stable 
period and highly orthogonal (86°)) to undertake this anal- 
ysis. We arbitrarily assign SPH particles with temperatures 
T<40000K to be 'cold gas', although our analysis is not 
sensitive to this choice. For the stellar disc, we concentrate 
upon stars we know to be rotationally supported which were 
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Figure 6. Time evolution of the alignment profile of the halo's 
minor axis with the simulation aligned to the polar disc. The 
black line is for the z-axis (parallel to the polar disc angular mo- 
mentum), the red line for the x-axis, and cyan the y-axis (ap- 
proximately parallel to the stellar disc angular momentum). The 
ordinate for each panel is the arccosine of the modulus of the 
projection of the minor axis along each of the principle axes. The 
i/-axis is defined as the angle of the minor axis relative to the 
z-axis of the angular momentum of the cold gas. 



formed between 6—6.5 Gyrs after the start of the simulation, 
in order to minimise contamination from the old centrally- 
concentrated disc+bulge stars, the stellar halo, and the polar 
disc itself. Using the circular velocity (v C i rc ) and the galacto- 
centric radius (r), we calculate the mass within that radius 
(M) as: 



M{< r) = 



G 



(3) 



where G is the usual gravitational constant. 

Using the circular velocity profiles of the two orthogonal 
discs we can probe the dark matter halo using observation- 
ally measurable galaxy properties. We show the velocity and 
enclosed mass profiles of the discs in Fig 7 (top and bottom 
panels, respectively). If there is more mass within a distance 
r of the centre of the galaxy then the halo is compressed in 
that direction. From the bottom panel of Fig 7, we can see 
that there is more mass along the stellar disc than the polar 
disc, which corresponds to a halo compressed orthogonally 
to the polar disc. At 3 kpc the difference in mass is ~70%, 
with the mass along the stellar disc being greater. Taking 
the cube root of the mass, to crudely estimate the axial 
ratios, gives 0.9, which matches remarkably well with the 
direct computation of the halo shape, which was found to 
be 0.92 (Fig. 14.11) . confirming the validity of using the two 
disc s to probe th e shap e of the halo. 

iBrook et al l (|2008l ) note a very high V/a relation for 
the stellar disc of 4.4 and 1.4 for the polar disc. This 
matches well w i th th e high V/a seen in observations, e.g. 
ISchweizer et all (ll983Ti who ca lculate a V/a of 2.2 for A0136- 
0801 and lSpavone et alj l|201ll ) who give V/a = 1.8, 2.15 and 
1.28 for UGC 7576, UGC 9796 and NGC 4650A respectively. 




Figure 7. Top panel: Line-of-sight velocity profile for the inter- 
mediate age stars (red line) and the polar disc cold gas (black 
line) at redshift z=0.17 (age=11.6 Gyr). Bottom panel: Enclosed 
mass as a function of radius based upon eqn[3] 



5 FORMATION AND EVOLUTION 
5.1 Angular Momentum 

In order to understand the processes that give rise to the po- 
lar disc, we examine the evolution of the angular momentum 
of the gas, stars, and dark matter. Although not directly 
related to the disc structures, the effect of dark matter is 
important because of the effect of tidal torques. 

In order to explore the formation processes of polar disc 
systems, we next traced the origin of the gas t hat lies within 
the sy stem and formed stars in the discs. iMaccio et al.l 
(2006) point out that unless the infall of gas is along the 
polar axis of the potential, it will be dragged to align with 
the stellar disc potential. When gas collapses under grav- 
ity it retains angular momentum and thus, the gas which 
comprises the polar disc and the stars of the stellar disc at 
2=0.17 should have a different origin. Two possible origins 
of the polar disc galaxy are suggested: first, the gas that 
comprises the two discs falls along two different trajecto- 
ries with two different angular momenta, or second, the last 
major merger was inclined to the galaxy and so re-oriented 
the existing gas and stars. New gas would then fall in along 
the old filament, but the early infall baryons would not be 
re-oriented. 

We align the simulation such that the polar disc at 
redshift 2=0.17 is aligned to the z-axis and use that ori- 
entation at each earlier timestep. The evolution of the an- 
gular momentum of the inner and outer dark matter halo, 
the gas, and the stars, is shown in Fig 8. This plot shows 
that the alignment of the angular momenta of the old stars, 
intermediate-age stars, gas, and dark matter experience a 
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Figure 8. Normalised angular momentum projected on the three 
principle axes of the simulation volume : black shows the z-axis, 
red the x-axis, and blue the y-axis. The top panel, labelled 'halo', 
is the angular momentum of the entire halo; the panel labelled 
'inner' is the angular momentum of the inner kpc of the halo; the 
panel labelled 'stars' is the angular momentum evolution of old 
(t<3 Gyrs) stars; the panel labelled 'cold' shows the trend for the 
cold gas. The grey area is the time of the last major merger. 

constant evolution with time after 6 Gyr. The polar disc 
galaxy evolves continuously into its z—0.17 (11.6 Gyr) posi- 
tion with no sharp transitions. 

The angular momentum of the dark matter halo (panel 
a) is established early and remains essentially constant for 
the rest of the simulation, although some tumbling can be 
seen in the z-axis. The angular momentum of the entire 
halo remains essentially the same throughout the history of 
the simulation, especially the x-component. This matches 
with the shape of the dark halo in the outer regions, which 
remains essentially constant. 

We see that the angular momentum of the inner region 
of the dark halo (panel b) and the stars are locked together 
because in this region the stars dominate the gravitational 
potential of the halo. 

We have previously discussed that the shape of the polar 
disc halo, and its orientation, remains fairly constant with 
time. We can use the axes of the best-fit ellipse at a radius 
of 0.25R V ir to calculate the cross-correlation of the gas and 
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Figure 9. The angle between the angular momentum of the cold 
gas (solid line) and the old stars (dashed line) and the major (top 
panel), intermediate (middle panel) and minor (bottom panel) 
axes of the dark matter halo. The shape of the halo is measured 
at 0.25R v i r . The dotted green lines show the time of the last 
major merger for reference and the red lines indicate the angular 
momentum which the galaxy would presumably have followed 
without the last major merger occurring. 



stars with the shape of the halo potential, as shown in Fig. 

i5~n 

The angular momentum of the polar disc aligns to the 
minor axis of the dark halo, as previously discussed, and the 
old stars align with the major axis within 2 Gyrs of the last 
major merger. The cold gas takes up to 4 Gyrs to move into 
close alignment with the minor axis. This may be suggestive 
of interactions between the dark halo potential which follows 
the dark halo shape and the galaxy gas and stars 

The old stars (Fig. I5.ll panel c) move orthogonally to 
the cold gas (Fig. I5.ll panel d) after the last major merger. 
During the last major merger there is a peak in the z-axis 
angular momentum of the stars which appears to have no 
long-term effect. This can be seen in both Fig. I5.1l and Fig. 
15. 1 1 At this point, the other two axes begin to twist into their 
final redshift 2=0.17 (age=11.6 Gyr) position. 

A striking feature of Fig. 15. H is that there is no difference 
in the orientation of the gas and stars along the intermedi- 
ate axis in any given time period (middle panel). The cold 
gas angular momentum becomes increasingly aligned to the 
halo minor axis (bottom panel). The most noticeable fea- 
ture from this plot is that the gas sees a temporary change 
in direction relative to the major axis (top panel) after the 
major merger, but recovers. The red line shows the assumed 
evolution of the gas angular momentum had the last major 
merger never occurred although this has not been tested. 
The old stars, however, do not recover or regain their old 
angular momentum. 

The angular momentum of the galaxy at early times 
is closely aligned with the intermediate axis of the halo in- 
ertial tensor but becomes increasingly divergent. The gas 
and stars are originally closely aligned but diverge after the 
last major merger. We see that initially the galaxy angu- 
lar momentum is orthogonal to the major and minor axes 
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and aligned to the intermediate axis. After the last major 
merger there is a change in the angle of both components 
relative to the major axis that the old stars maintain. The 
infalling gas returns to orthogonality within several Gyr. 
Without the major merger it seems likely that the gas and 
stars would remain orthogonal to the major axis. There is, 
however, an evolution in the angular momentum relative to 
the intermediate axis, such that the gas and stars angular 
momentum is initially along this axis but then increasingly 
diverges. Had the major merger not disrupted the system it 
seems likely that there would be a continuous evolution of 
both components until the galaxy a ngular momentum aligns 
to the minor axis in agreement with lBailin et all (|2005h This 
can be assumed from Fig. 15.11 The recovery time suggested 
corresponds to a feature of the total angular momentum of 
the cold gas. Although not shown in Fig. 15. li the last major 
merger appears to result in an increase in the total angu- 
lar momentum of the cold gas and if followed by a trough 
which lasts from 4 Gyr to 6 Gyr. After 6 Gyr the angular 
momentum of the gas rapidly increases again. The angular 
momentum of the old stars does not change in magnitude 
during this collision. 



5.2 Formation Process 

The formation process of the galaxy is intimately connected 
with the evolution of the angular momentum described in 
the previous section. This is because the changing alignment 
of the galaxy with respect to the large scale environment is 
what gives rise to the polar disc properties. 

Another manifestation of misaligned gas and stellar 
discs can be seen in the phenomenon of warps. It is perhaps 
reasonable to assume that warps and polar discs might have 
a common origin. In this section we first briefly examine the 
possibility that a polar disc is an extreme type of warp. We 
then t est the formation scenario presented in I Brook et al.l 
(2008) by tracing the angular momentum evolution of the 
infalling gas, and the direction of the last major merger. 

Warps are now thought to be due to the interac- 
tion between infal ling gas and misaligned hot gas (e.g., 
iRoskar et all l|201Ch ). In our simulated system, though, the 
hot gas is aligned to the polar disc. Warps begin at the trun- 
cation of the stellar disc and the g as is only misaligned to 
the s tellar disc outside this radius |Sanchez-Blazquez et al.l 
120091 ). The polar disc at z=0.17 (11.6 Gyr) is out of align- 
ment at all radii and (essentially) perfectly orthogonal, sug- 
gesting that polar discs and warps are distinct processes. At 
first glance, there seems to be no significant similarity be- 
tween warps and polar discs. However, there is evolution of 
the direction of the polar disc galaxy, as a whole, after the 
formation of the polar disc (Fig. I5.1|l , possibly does to tidal 
torques. 

We identified gas that comprises the polar disc as all 
gas that is cooler than T=40000K within the inner 10 kpc 
of the host halo at z=0.17 when the polar disc is most promi- 
nent. This lower temperature threshold does not change the 
direction of the angular momentum of the gas but results 
in a thinner, better defined, disc structure. We track these 
particles to their origin, then repeat the process tracing the 
gas that formed stars prior to the last major merger at 2=2. 
When we plot the positions of these sets of particles there is 



little obvious difference between the two, which one might 
expect if the origin of the gas was different filaments. 

If we examine an output from before the last major 
merger, we see that there is a broad range of angular mo- 
mentum in the infalling gas which becomes systematically 
more aligned to the z-axis with time. The angular momen- 
tum of the infalling gas is surprisingly coherent with dis- 
tance from the galaxy. The total angular momentum of the 
gas which forms the polar disc stars is tightly aligned with 
the z-axis from ~4 Gyrs onwards. It is a similar story if we 
exclude the gas which is already in the inner region (within 
the inner 1 kpc of t he galaxy). This appears to confirm the 
iMaccio et alj (|2006l ) hypothesis regarding the importance of 
cold flows along filaments for polar disc galaxy formation. 

Another feature of the polar disc structure is that dur- 
ing intermediate times the star formation in the old stellar 
disc is fuelled by newly infalling gas. The reason for this is 
that when gas falls into the innermost region of the galaxy 
it is strongly affected by the already formed stars, which 
changes the angular momentum of the gas and brings it 
into alignment with the stellar disc. In the inner region the 
potential is dom inated by the stars. An obser vational confir- 
mation of this is lCombes fc Arnaboldil (Il996l ) who find that 
no dark halo is needed to reproduce the rotation curve of 
the stellar disc in NGC4650A. 

This behaviour no longer occurs after about 9 Gyr. As 
the entire galaxy is gradually rotating to bring it into align- 
ment with the galaxy halo shape; it may be that the gas is at 
a steep enough angle to the stellar disc so that it is no longer 
dragged into alignm ent with the stellar disc at later times. 
This is suggested bv lMaccio et al.l (|2006T ) as to why infalling 
gas needs to be approximately orthogonal for a polar disc to 
be stable. 

We also see that a considerable fraction of the gas which 
forms the polar disc is aligned to the z axis when it is still 
far from the galaxy, but by no means all the gas is tightly 
aligned, and interactions within the galaxy halo are impor- 
tant in finally aligning the gas. 

The scenario favoured bv lBrook et alj (|2008l ) suggested 
that the polar disc galaxy is formed if the last major merger 
has a total angular momentum orthogonal to the infalling 
gas. This interaction re-orientates the old gas and stars while 
the angular momentum of newly infalling gas remains un- 
changed. The coincidence of the last major merger with the 
more unusual behaviour of the angular momentum (Fig 8) 
of the galaxy is suggestive that it is essential to forming the 
galaxy. The mass of stars in the polar disc galaxy at z—2 is 
3.2x 10 9 M Q and the merging satellite galaxy is 1.5x 10 9 M Q , 
making the interaction at least 1:2. 

The angular momentum of the collision is calculated 
by assuming that each galaxy is a point sitting within the 
reference frame of the infalling gas. The relative tangential 
motion of the collision is anti-clockwise. A schematic exam- 
ple of the galaxy-galaxy interaction is shown in Fig 10. By 
examining the internal angular momentum of the interact- 
ing galaxies we found that most of the angular momentum 
of the interaction is stored in the bulk motion of the galaxies 
rather than in internal motions. The angular momentum of 
the satellite's stars is aligned to the angular momentum of 
the merger. 

The angular momentum of the collision is along the x- 
axis in the chosen reference frame. This is orthogonal to the 
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Figure 10. Projected schematic of the last major merger at red- 
shift 2=2 (age=3.3Gyr), just before the stars merge. The black 
star is the central galaxy, while the red star is the galaxy which 
merges. The blue and green lines are the projected velocity vec- 
tors. Distances are in kpcs. 



angular momentum of the gas falling into the halo from the 
filament which is aligned to fall in along the z-axis. As the 
angular momentum of the system must be conserved, we 
suggest that this provides the initial kick which re-orients 
the old stars away from the initial orientation, see Fig. 15.11 
This figure suggests that the major merger alters the be- 
haviour of the galaxy away from a matched gas and stellar 
disc which reorientates from being aligned to the halo inter- 
mediate axis to the minor axis. 

Much of the gas that forms stars between 4 and 6.5 Gyr 
has its origin in the satellite galaxy involved in the last major 
merger, and thus contributes greatly to the alignment of the 
stellar disc. Specifically, 40% of the stars formed between 4 
and 6.5 Gyr were produced from gas which fell in with the 
last major merger. 

The motion of the system after the major merger may 
be due to slow tumbling or torquing by the halo the halo as 
mass continues to be accreted. This is important because the 
stellar disc eventually aligns with the y-axis of the system. 



5.3 Metallicity and Age 

We derived the radial (mass-weighted) metallicity gradients 
of the discs associated with our simulated polar system by 
placing narrow slits along each component, ignoring the in- 
nermost 1 kpc of both, to avoid bulge contamination and 
the intersection point between the two. The overall stellar 
metallicity of the stellar disk is roughly twice (four times) 
that of the polar stellar (cold gas) disc. 

Fig [TT] shows the radial metallicity profiles for the 
intermediate-age stars of the stellar disk (black) and the 
cold gas of the polar disc (red). The metallicity gradi- 
ent of the latter is very flat (~— 0.005— > — 0.01 dex/kpc) 
consistent with the recent obs e rvatio n of flat gradients in 
NGC 4650A, bv lSpavone et all (|2010T ). The metallicity gra- 
dients of the stellar and polar discs are —0.09 dex/kpc and 
—0.04 dex/kpc, respectively. This matches the results of 
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Figure 11. Radial mass- weighted metallicity profile for the stars 
formed in the primary stellar disc (red), and the cold gas of the 
secondary polar disc (black), along 1 kpc thick slits. The stel- 
lar disc truncates at ~4 kpc; beyond that point, star particles 
are kinematically associated with the stellar halo. The polar disc 
truncates at ~8 kpc. The dotted lines are the best fit lines to 
the metallicity profiles, where [Z/Z ] = m * r + c and m is the 
gradient in dex/kpc, and c is the y intercept. 



ISpavone et al] l|201lf ) who find that the profile for UGC 7576 
is also extremely shallow. 

The negative metallicity gradient seen in Fig II II is char- 
acter is tic of inside-out disc g ro wth |Matteucci fc Francois! 
1 19891: iMatteucci et al] 1 19891 : iBoissier fc Prantzosl Il999l : 
iPilkington et al.l 12012). consistent with t he fo rmation sce- 
narios suggested by iMaccio et al. I (|2006T ) and iBrook et al.l 
(2008) for their simulated polar disc systems. We did how- 
ever note that the gas in the polar disc initially forms a ring 
at larger radii and then later 'fills in' the inner region. This 
is because while the polar disc is forming, the gas in the 
inner region is inclined to the old stellar disc rather than 
the developing polar disc, depleting the inner region of the 
polar disc structure. We suggest that the abundance gradi- 
ent is indicative of the stellar disc enriching the gas in the 
inner region of the polar disc rather than polar disc stars 
alone. The lower metallicity polar disc is consistent with its 
younger age sampling an earlier part of the age-metallicity 
relation. It is consistent with a polar disc being preferen- 
tially 'polluted' by the infall of 'primordial' gas flowing into 
the system at late times along filaments. When we track the 
origin of this gas below, we confirm that this is indeed the 
case. 

In terms of global metalli city, this polar disc g a laxy is 
similar to those surveyed by ISpavone et al.l l|2010l . 1201 if ) . 
Specifically, Z/Zq — 0.2 for our galaxy, matching the metal- 
licity of NGC4650A, higher than UGC 9796 (Z/Z =O.l ) 
and lower than UGC 7576 (Z/Z =O.4). 

The radial age profile of the stars along each disk shows 
that the stars are significantly younger in the polar disc than 
in the stellar disc (Fig I12|l . The older spheroid component 
tends to dominate at both low and high radii, driving the 
declines seen at the extrema of both curves. If we remove 
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Figure 12. The stellar age profile along the polar (black) and 
stellar (red) discs. i/ orm is the mean formation time of stars 
at a given radius. The stellar (polar) disc truncates at ~4 kpc 
(~8 kpc). The black dotted line shows the star formation time 
with radius for the polar disc, discounting old stars (tf orm < 
AGyr). 



stars that formed in the first ~4 Gyrs of the simulation, 
we find a flat stellar age gradient in both the stellar and 
polar discs. The mean formation time of stars in the stellar 
(polar) disc are 5—6 Gyrs (~10 Gyrs) from the start of the 
simulation. 

Loo king at Fig. EH . we se e a much higher SFR than is 
found bv lSpavone et all (|201l[ ). However, their SFR matches 
that for our polar disc at 9 Gyr or about 1 Gyr after it 
begins to form. This is a similar age to the estimated age of 
the polar disc in NGC 4650A. 

As part of our examination of the temporal evolution 
of the angular momentum, we looked at an alternative ap- 
proach to decomposing the two discs, aligning the simulation 
volume using the cold gas distribution and calculating the 
angular momentum of each star particle. We can then plot 
the ratio of the z- and ^/-components of the angular momen- 
tum to the total angular momentum, as shown in Fig 1131 It 
is noticeable that the angular momentum of stars in the po- 
lar disc is tightly correlated with the gas, such that all stars 
with a Jz/Jt greater than 0.9 are considered part of the 
polar disc. For the stellar disc, the distribution is broader, 
suggesting stars with J y /Jt greater than 0.75 are part of the 
stellar disc. 

The above provides an age- independent method to iden- 
tify the two discs using angular momenta without the age 
cuts used previously. This allows us to see how the discs 
grow, we bin the stars according to their distances along 
their respective discs. If we repeat the calculation of the stel- 
lar age profile using this new method we can see that the 
slope at radius < 4 kpc and > 6 kpc is apparently due to 
contamination from the bulge and stellar halo respectively. 
This shows that there is no significant difference in the ages 
of stars at different radii, despite the presence of metallicity 
gradients. 

Fig. 15.31 panel (a), shows that the star formation his- 
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Figure 13. The angular momentum profiles of stars along the x 
(red), y (blue), and z (black) axes, where the cold gas has been 
aligned with the z-axis. 



tory of the polar disc changes slightly with radius, while 
panel (b) show the same plot for the older disc. In the po- 
lar disc, panel (a) shows that there is a population of stars 
which forms early, contaminating the sample. Indeed, if we 
remove those early forming stars (tf orm < 4Gyr) from the 
age profile then the age gradient is almost completely flat, 
as if we were using the decomposition method discussed di- 
rectly above. The black dotted line in Fig 12 shows this 
trend. In Fig. 1111 we show the gas metallicity gradient for 
the polar disc. When we trace the gas back in time, and 
identify it just before it enters the inner region (radius < 12 
kpc) we find that the metallicity gradient of the infalling gas 
is considerably flatter than at z=0.17. Additionally, we find 
no particular gradient in the time at which the gas entered 
the inner region. This suggests that the reason for the metal- 
licity gradient at z=0.17, is the enrichment of the gas in the 
inner region due to a higher star formation rate at small 
radii, which is uniform across time. The metallicity gradient 
is, apparently, not indicative of an inside-out formation of 
the polar disc. 

The metallicity distribution of the polar disc stars has 
the same gradient as the cold gas if we remove the very old 
stars, but the stars are 0.15 dex less rich at any given radius. 
If we include the old stars then we find that the inner region 
is considerably more metal rich. 

We conclude that the lBrook et al.l (|2008l ) polar disc is 
comparable to the p o lar disc gala x ies ob served in nature 
ISpavone et all i|2010l ); llodice et al.1 (|2002T l. The polar disc 
structure is older than the value derived for NGC 4650A by 
Spavone et al., but has a comparable mean metallicity. The 
stellar disc is more metal-rich than the younger polar disc, 
and the metallicity gradient for the polar disc is comparably 
shallow to that observed in NGC 4650A. The polar disc 
does not produce stars for ~3— 4 Gyr after the last major 
merger, consist ent, for example , with the picture proposed 
for UGC 9796 (|Cox et aij|2006h . 
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(a) Polar disc (b) Stellar disc 



Figure 14. Age distribution of stars at different radii along each disc. The lines are colour-coded as stated in the inset to each panel. 
| (a) | shows the polar disc stars, while [(b)] shows the stellar disc stars 



6 CONCLUSIONS 

We have extended the prelimi nary analysis of the polar disc 
structure first introduced by I Brook et al.l f200&) . We have 
examined the properties of the dark matter halo shape, the 
formation process of this object and the metallicity and stel- 
lar age gradients of the discs. 

Ilodice et al.l £2003) point out that polar disc galaxies 
provide an unparalleled chance to study the shape of the 
inner region dark matter haloes using observations in a way 
normally only possible in simulations. The two orthogonal 
discs have a circular velocity defined by the potential, which 
is dominated by dark matter. By comparing the circular 
velocity of orthogonal discs we estimate the flattening of 
the dark matter in the plane where both discs are seen edge 
on and find that the halo is flattened in the direction of the 
polar disc, as found by Iodice et al. Specifically the axial 
ratio is found to be c/a=0.9. Directly measuring the dark 
matter halo in the simulation we find the same result and 
compare find the axial ratios are extremely similar using 
both methods, i.e. c/a=0.93. The value quoted is for the 
inner region of the halo, at the radius of the polar disc itself. 
In the outer region of the disc the sphericity declines to 0.67. 
We concentrated on the inner region in this work because 
that is the region probed by the baryons. 

This value is closely matched by simulations of dark 
matter haloes, su ch as the galaxy size h alos with cool- 
ing simulated by iKazantzidis et all |2004) and shown in 
their figure 2, which shows an almost spheroidal inner re- 
gion, with a sphericity of 0.9 at 10% if the virial radius, 
and an outer r egion of c/ a=0.6, confirm i ng th e work of 
ISpringel et all d2004l) and lAllgood et all (|2006l ) amongst 
others. lAllgood et all I^OOa T has categorised halo shapes 
from cosmological simulations using statistical samples. Us- 
ing their results we expect our halo to have a c/a ratio of 
approximately 0.65 at z=0.17 which is consistent with the 
polar disc galaxy halo; the value of c/a for the 'classical' disc 
shown in Fig. 14.11 at the viral radius is 0.7. The shape pro- 



file, both in angle and sphericity remain stable from outside 
the inne r region to the v irial rad ius. This differs from the 
work of iBailin et all (|2005l ) and lAllgood et ail (|2006h who 
find a change of angle and shape with radius. This feature 
is unusual but po lar discs are unusua l objects and this may 
be characteristic. iMaccio et ail l|200rj) did not see the same 
level of coherence in their polar disc galaxy, however. 

Compared to observational estimates, e.g. [Sackctt 

(1999), however, these halos appear very spherica l. Stud - 
ies of polar ring/dis c gala x ies (see Sackett et"al1 (1 19941 ). 
ICombes fc Arnaboldil (|l996l V Ilodice et al l ll2003l)~ sugg est a 
highly flattened halo of between 0.3 and 0.4. ISackettl (| 19991 ) 
does, however, provide some suggestions as to why this dif- 
ference occurs. Estimates of c/a for halos based on different 
observations vary wildly depending on the chosen metric, 
such that c/a — 0.5 ± 0.2. In the outer region our halo falls 
within this range, 0.65, while in the inner region (R < 15 
kpc) where the observational estimates probe we see c/a=0.9 
to 0.75. The simulatio n outcome is also close to the result of 
iHoekstra et~aH (|2004l ) based on weak lensing who estimate a 
halo sphericity of 0.67lg q?i the upper range of which is close 
to the sphericity of our halo at a radius of 14 kpc. As for the 
extremely elliptical shape found for polar ring/disc galaxies, 
the methods used m ay simply not give good estimates, in- 
deed, ISackettJ (|l999l ) notes the difficulty in measuring halo 
sphericity observationally. 

We have concluded that the polar disc of this system 
is not caused directly by a merger or stripping. There is 
no single merger or stripped satellite which contains enough 
gas to form the massive polar disc. By process of elimina- 
tion, the pola r disc must form from inflowing cold gas, as in 
M accio et ah! l|2006l ). We attempted to identify the mecha- 
nism by which this cold gas flows into the system and gives 
rise to the polar structure, asking, why the filaments should 
change direction relative to the old stellar disc. The struc- 
ture of the polar disc seems to be a result of the direction of 
gas infall from the filament. We have shown that there is not 
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a single moment when the polar disc 'forms' but, instead, 
there is a constant evolution of the angular momentum of 
the infalling gas and stars. 

The angular momentum of the last major merger is al- 
most exactly orthogonal to the filament and we feel that 
this is the most likely origin of the polar disc st ructure. In 
this scenario, put forwards bv lBrook et al.l |2008), the major 
merger reorientates the angular momentum of the old stars 
and the cold gas already present in the disc. Gas continues 
to flow in from the large scale structure unaffected by this 
merger and so builds up the polar disc. 

We see that the major merger affects the stars most 
strongly without any significant influence on the dark halo, 
and the cold gas recovers its old angular momentum within 
a Gyr. However, this is not the whole story as subsequent to 
the major merger there is further evolution in the orientation 
of the entire galaxy system, potentially due to tumbling or 
torques from the halo. The growing alignment of the polar 
disc with the intermediate axis of the dark halo does not 
seem affected by the major merger. 

We also provide a basic comparison of the polar disc 
galaxy to a classical disc galaxy to better isolate the origin 
of the extreme behaviour of the polar disc. We also confirm 
that the polar disc exhibits the characteristics of a disc built 
up from the inside out, both from the metallicity gradient 
and the stellar age profile. Clearly though, any conclusions 
for their formation, when based upon a single simulated sys- 
tem, must be taken as tentative, at best. We need to (and 
will) repeat this analysis on a statistically significant sample 
of simulated polar disc systems. 
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